Astronomy & Astrophysics manuscript no. aal8890 
April 3, 2012 



©ESO2012 



Letter to the Editor 



Damping rates of solar-like oscillations across the HR diagram 

Theoretical calculations confronted to CoRoT and Kepler observations 

K. Belkacem 1 , M.A. Dupret 2 , F. Baudin 3 , T. Appourchaux 3 , J. P. Marques 41 , and R. Samadi 1 



(N 

o 

(N 

< 

CO 

6 

S-H 
-t— > 

> 

m 
i> 

en 

o 



X 



1 LESIA, UMR8109, Universite Pierre et Marie Curie, Universite Denis Diderot, Obs. de Paris, 92195 Meudon Cedex, France 

2 Institut dAstrophysique et de Geophysique, Universite de Liege, Allee du 6 Aout 17-B 4000 Liege, Belgium 

3 Institut d'Astrophysique Spatiale, CNRS, Universite Paris XI, 91405 Orsay Cedex, France 

4 Georg-August-Universitat Gottingen, Institut fur Astrophysik, Friedrich-Hund-Platz 1, D-37077 Gottingen, Germany 



April 3, 2012 



ABSTRACT 



The space-borne missions CoRoT and Kepler are providing a rich harvest of high-quality constraints on solar-like pulsators. Among 
the seismic parameters, mode damping rates remains poorly understood and are thus barely used to infer the physical properties 
of stars. Nevertheless, thanks to the CoRoT and Kepler spacecrafts it is now possible to measure damping rates for hundreds of 
main-sequence and thousands of red-giant stars with unprecedented precision. By using a non-adiabatic pulsation code including a 
time-dependent convection treatment, we compute damping rates for stellar models that are representative of solar-like pulsators from 
the main-sequence to the red-giant phase. This allows us to reproduce the observations of both CoRoT and Kepler, which validates 
our modeling of mode damping rates and thus the underlying physical mechanisms included in the modeling. By considering the 
perturbations of turbulent pressure and entropy (including the perturbation of the dissipation rate of turbulent energy into heat) by the 
oscillation in our computation, we succeed in reproducing the observed relation between damping rates and effective temperature. 
Moreover, we discuss the physical reasons for mode damping rates to scale with effective temperature, as observationally exhibited. 
Finally, this opens the way for the use of mode damping rates to probe turbulent convection in solar-like stars. 
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1. Introduction 

The space missions CoRoT and Kepler are providing accu- 
rate observations of solar-like oscillations of hundreds of main- 
sequence stars and thousands of red-giant stars. Therefore, we 
now have access to seismic properties, such as mode frequen- 
cies, amplitudes, and linewidths, for a large and homogeneous 
sample of stars. Among them, the physical mechanisms gov- 
erning mode linewidths (T), which is linearly related to mode 
damping rates {rf) such as F — rj/n, are still poorly understood. 
Many studie s have attempted to model the damping rates of solar 
mode s (e.g.. lBalmforthlll992t iDupret et alJl2006al IXiong et alJ 
l2000l) . All those modelings have had difficulties in reproduc- 
ing the solar observations except when using free parameters, 
and the under lying physical mechanisms remain unclear (e.g., 
lHoudekll2008l) . 

The scaling relation of mode linewidths for other stars 
thus provides an important additional constraint on the mod- 
eling. A usual way to investigate mode parameters adopts 
scaling relations bet ween asteroseismic quanti t ies and stel- 
lar parameters (e.g. Kieldsen & Beddind 1995 : Samadi etal. 



2007tlMosser et al.i20irilBaudin et alj20l"la b; Belk acem et al 



20111) . For mode linewidth s, scaling relat i ons ha ve been inves- 
tigated only very recently. IChaplin et alJ (120091) used ground- 
based observations and proposed that mode linewidths follow 
a power-law dependence on effective temperature, i.e. F oc T 



c-ir 



where F is the mode linewidth and T e g the effective temperature. 



Howev er, those conclusions were challenged by iBaudin et al.1 
(l2011allbh using a homogeneous sample of stars observed by 
CoRoT. They found that a unique power-law cannot describe the 
entire range of effective temperature covered by main-sequence 
and red-giant stars and concluded that mode linewidths of main- 
sequence stars follow a power-law of T 1 ^* 2 , while red-giant stars 
only slightly depend on effective temperature (7!^ 9 ). This re- 
sult, for main-sequence and sub-giant stars, was later extended 
by Kepler observations (I Appourchaux et aTll2012l) . 
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From a theoretical point of view, IChaplin et al.l d2009l) pre- 
dicted a power-law of T oc r 4 ff , which d isagrees with CoRoT 
and Kepler observations. iHoudekl (1201 2h attributed the failure 
of this theory to a missing physical mechanism and proposed 
mode scattering as a possible solution. This disagreement em- 
phasizes that our understanding of mode linewidth is still in 
its infancy and that a first necessary breakthrough would be 
to reproduce the CoRoT and Kepler observations. Hence, in 
this letter we adopt the f ormalism of iGrigahcene et alJ (120051) 
and IDupret et al.1 d2006bl) . This allows us to confront theoreti- 
cal computations with observations from both the CoRoT and 
Kepler spacecrafts, and consequently to understand the strong 
dependence of mode linewidth on effective temperature and val- 
idate the theoretical modeling. 
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2. Computation of theoretical damping rates 

2.1. The grid of stellar models 

We used a grid of stellar models for masses between M = 
lAf Q and M - 1.4 M Q from the ZAMS to the tip of 
the red giant branch, which are typical of observed solar- 
like pulsators. The grid was obtained using the stellar evo- 
lution code CESAM2k dMorel & Lebretonl 12008k The atmo- 
sphere was computed assuming a grey Eddington approxima- 
tion. Convection was included according to the Bohm-Vitense 
mixing-length (MLT) formalism, with a mixing-length pa- 
rameter o- = 1.6. T he initial chemical composition follows 
lAsplund et alj d2005l) . with an helium mass fraction of 0.2 485. 
We used the OPAL equation of s tate (R ogers et alj|1996l) and 
opacities (Ialesias & Rogers 1996), complemented, at T < 10 4 
K, by the Alexander & Ferg uson ! d 19941) opacities . We u sed the 
NACRE nuclear reaction rates from lAngulo et alj d 1999b except 
fo r the 14 N + p reactio n, where we used the reaction rates given 
in llmbriani et alj ( l2004l) . 

2.2. The non-adiabatic oscillation code 

Damping rates were computed using the non-adiabatic pulsation 
code MAD (Dupret 2001), which includes th e time-dependent 
convec tion (TDC) treatment described in iGrigahcene et alj 
d2005l) . This approach takes into account the role played by the 
variations in the convective flux, the turbulent pressure, and the 
dissipation rate of turbulent kinetic energy. This TDC a pproach 
is a no n-local formulation of c onvection based on th e lGabriell 
d!996l) formalism explained in iDupret et all (l2006blfah . In this 
framework, non-local parameters related to the co nvective flux 
and the turbulent pressure are chosen such as in IDupret et al.l 
d2006bl) so that it fits the solar 3D numerical simulation. 

In addition, it involves a parameter p 1 , which takes com- 
plex values and enters the closure term of the perturbed en- 
ergy equation. This parameter was introduced to prevent the oc- 
currence of non-physical spatial oscillations in the eigenfunc- 
tions (see IGrigahcene et al.ll2.Q05l . for details). To constrain this 
parameter, we adopt the following strategy: (3 is adjusted so 
that the frequency of the depression of the damping rates (see 
Fig. [TJ coincides with v max , where v max is computed by using 
the linear relation between v max and the cut-off frequency (e.g., 
iKjeldsen & Beddinglll995l) . We note that the bond between v max 
and the frequency of bottom of the secon d ascending branch of 
the damping rates is observed for the Sun dBelkacem et alj201 ll) 
and m ore generally for solar-like pulsators dAppourchaux et all 
12011 . 



2.3. Dominant contributions of the damping rates 

The integral expression of the damping rates can be written as 
dGrigahcene et alj|2005l) 
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where a> is the mode frequency, / the mode inertia, (Ft, - 1) = 
(d In To/d lnpo)s, Sp the Lagrangian perturbation of density, 5S 
is the Lagrangian perturbations of entropy, 8P tur b the perturba- 
tion of turbulent pressure, To the mean temperature, po the mean 
density, and the star denotes the complex conjugate. 

The first term of Eq. (Q~|) includes the contributions of the per- 
turbations of the radiative and convective fluxes, as well as the 
perturbation of the dissipation rate of turbulent kinetic energy 



1.0000 1 



0.1000 r 



3. 0.0100 r 



0.0010 r 



0.0001 




2.0 2.5 3.0 
v [mHz] 



1.0 



Fig. 1. Damping rates versus mode frequency, computed for a 
model of a mass of M - 1.1 M and an effective temperature 
of r e ff = 6026K on the main sequence. The damping has been 
computed as described in Sect. 12.21 The solid line corresponds 
to the total damping rate, while the dash-dotted line corresponds 
to the contribution of the perturbation of turbulent pressure and 
the dashed line to the absolute value of the contribution of the 
entropy perturbations (see Eq.[TJ. The vertical dotted line corre- 
sponds to v max . 



into heat. We note that this term corresponds to the non-adiabatic 
part of the gas pressure perturbation, which is found to be nega- 
tive (thus a driving contribution). The second term of Eq. (Q]i rep- 
resents the perturbation of turbulent pressure. Figure [TJ displays 
the mode damping versus the mode frequency as well as the two 
contributions expressed in Eq. ([TJ. It turns out that the contribu- 
tion of turbulent pressure dominates the damping and is partly 
compensated by the contribution of entropy. Both contributions 
have roughly the same order of magnitude, hence the total mode 
damping is small compared to the absolute values of both the en- 
tropy and turbulent pressure contributions. In addition, we note 
that the depression (or plateau) of the damping rates is the result 
of the maximum compensation between the two contributions. 
This is the case for all the models, from the main-sequence to 
the red-giant phases. 



3. Comparison with CoRoT and Kepler observations 

3.1. Data set 

We considered CoRoT and Kepler stars for which linewidths 
had been accurately measu red. We first c onsider ed the CoRoT 
observations described in iBaudin et alj d2011alfbh . The mea- 
sured linewidths for main-sequence stars come from the 
CoRoT seismological field. More pre cisely, we used HD 49933 
dBenomar et all I2009I). HD1814 20 dBarban et alj l2009l) . and 
HD49385 dDeheuvels etai]|2010h . We also used the results fo r 
HD50890 (Baudin et al.ll2012l) . HD1 81907 dCarrier et all20ld) . 
and the Sun. For CoRoT red-giants, IBaudin et al.l d201 lallbl) fit- 
ted 235 red giants observed in the exo-field of CoRoT for about 
142 days. 

We also considered Kepler observatio ns of 42 stars ac- 
quired over nine months as presented by lAppourchaux et aD 
(|2012|) . Those stars were analyzed by several groups and mode 
linewidths at the maximum of height in the power spectrum were 
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Fig. 2. Mode linewidths (normalized by the solar value, r sun = 
0.95/jHz) versus effective temperature. The squared symbols 
represent theoretical calculations computed as explained in 
Sects. l2~Tl and l2. 21 The triangles cor respond to the observations 
of main-sequence stars derived by lAppourchaux et all d2012l) 
from the Kepler data (with their 3-cr error-bars). The dots corre- 
spond to the observations of red giants (with T e g < 5200 K) and 
main-sequence (w ith r e ff > 5200 K, w ith their 3-cr error-bars) 
stars as derived bv lBaudin et al.l (l2011allbl) from the CoRoT data. 



derived. We note that lAppourchaux et al.l(l2012l) considered only 
main-sequence stars and sub-giants but not red giants. 

3.2. Results 

A direct comparison between observations and theoretical com- 
putations is provided in Fig. [2] for both Kepler and CoRoT ob- 
servations. For main-sequence stars, there is an overall agree- 
ment between the theoretical computations and both CoRoT 
and Kepler observations within the observational error-bars. 
However, the observed Kepler linewidths are smaller than pre- 
dicted by theory and this is particularly the case for stars with 
high effective temperatures. While this result can be related to 
modeling deficiencies, there are also several possible observa- 
tional uncertainties. For instance, the determination of effective 
temperature is subject to i mportant uncertainties in t he estima- 
tion of the mode line width ( lAppourchaux et al1l2012l) . 

For red-giant stars, our theoretical computation is in overall 
agreement w i th the r esults derived for CoRoT observations by 
iBaudin et al.1 (1201 lallbh . We note however that possible uncer- 
tainties can also come from either the determination of effective 
temperatures or a bias due to the o bservation duration and thus 
the limited frequency resolution in IBaudin et al.l (l2011albh . An 
extended investigation of damping rates dedicated to red-giant 
stars would thus be desirable in a future work, and allow us to 
draw conclusions about the dependence of effective temperature 
to damping rates for red-giants. In addition, such a work will 
benefit from ongoing Kepler observations that will provide us 
red giants at low effective temperatures (T e ff < 4200 K). 

This overall agreement with both CoRoT and Kepler ob- 
servations demonstrates that the main physical picture is well- 
reproduced by the modeling. However, an extension of this com- 
parison using sub-giants and high luminosity red-giant stars ob- 



served by Kepler is desirable in the future to ensure that we have 
a homogeneous sample of stars across a large range of effective 
temperatures. 

4. Discussions 

We have discussed the strong dependence of mode damping 
rates on effective temperature and shown that it can be under- 
stood by simple arguments. Secondly, we have investigated the 
possible origin of a different scaling for the damping rates of 
main-sequence and red-giant stars. 

4.1. Mode linewidth versus effective temperature 

We start from Eq. (Q]) and note that there are obviously two im- 
portant quantities, the work integral and the mode inertia. The 
first depends on the phase difference between the mode com- 
pression (perturbation of density) and perturbation of pressure 
(gas and turbulent pressure). Hence, it is the non-adiabatic part 
of pressure fluctuations (6S /c v ) that mainly determines the mode 
damping. Accordin g to the mode energy equation (see, for in- 
stance, Eq. (A2) of iBelkacem et al1l201 lb . one can expect that 
the work integral dimensionally scales as the ratio L/M. This is 
confirmed by Fig. [3] (top), which shows that the work integral 
follows 



* 7oc Im) 



(2) 



In contrast, mode inertia (I) does not depend on mode energy 
leakage but on the star's static structure and more precisely on 
the properties of the uppermost layers. Hence, one can expect 
mode inertia to scale with the surface gravitjQ as in the case 
shown by Fig. [3] (bottom). More precisely, 



I CCg 



-2.4 



(3) 



Now using Eq. ((2j and Eq. (0 and further noting that L/M oc 
Ttglg, one easily finds that mode damping depends mainly on 
effective temperature, such that 



riKTT 



-0.3 



(4) 



Such a crude analysis is unable to reproduce the precise shape of 
the mode line-width with effective temperature. However, it al- 
lows us to explains qualitatively the strong dependence of mode 
damping rates on effective temperature. It turns out that the de- 
pendence on the effective temperature is a result of a compensa- 
tion between the work integral and mode inertia. 

4.2. Main-sequence and red-giants stars: two physical 
regimes? 



Several authors dAndoet al.l 120101: IBaudin et ail 1201 lallbh have 
shown by using both ground-based observations and photomet- 
ric observations from space that there is an apparent change in 
the behavior of the mode damping rates between main-sequence 
and red-giant stars. They therefore suggested that there is a 
switch between two different physical regimes. However, our 
theoretical results suggest that the physical mechanisms for 
mode damping remain the same from main-sequence to red- 
giant stars. 

However, the power-law of mode damping with effective 
temperature is very sensitive to the way v max is selected. This 



1 Note that mode inertia also scales with the dynamical timescale 
■\/(GM/R 3 ) with almost the same dispersion as for the surface gravity. 
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Fig. 3. Top: Product r\ x / (see Eq. ([TJ for details) as a function 
of the ratio L/M (where L and M are the luminosity and mass, 
respectively). Bottom: Mode inertia versus surface gravity g. 



is illustrated by Fig. [4] For instance, a bias in selecting v max can 
arise because in red giants some non-radial modes have almost 
the same widths as radial ones dDupret et alJl2009f) . For main- 
sequence stars, the mode density is important so that an error 
of, say, one radial order in the selection of v max has negligible 
effects. In contrast, for red giants, this shift can lead to the se- 
lection of a mode outside the depression of the damping rates. 
Hence, we conclude that the slope of mode linewidth, for red- 
giants, with effective temperature is sensitive to the way in which 
the frequency of the maximum height is de rived in each power 
spectr um. In turn, it could also explain why lAppourchaux et all 
(1201 2l) found different results from their fitting when the fre- 
quency of the maximum height or amplitude was selected. 



5. Conclusions 

In contrast to previous work, we have shown that the theo- 
retical computation of mode damping rates agree with both 
Kepler and CoRoT observations from main-sequence to red- 
giant stars. This result suggests that the main physical picture 
used to model mode damping successfully reproduces the obser- 
vations. Moreover, we have been able to understand the strong 
dependence of mode damping rates on effective temperature by 
pointing out that it is the result of a compensation between the 
work integral and mode inertia. 
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Fig. 4. As for Fig. [2] Double cross symbols correspond, as de- 
scribed in Sect. [3] to damping rates computed for the closest 
mode in frequency to the frequency of the maximum v max , minus 
one radial order. The crosses correspond to the same calculation 
but adding one radial order for the selection of the frequency of 
the maximum. This procedure is intended to mimic the effect of 
possible biases in the observational determination of v max . 



In contradiction again to previous results our theoretical 
computation of mode damping rates with the same physics de- 
scription are in agreement with both Kepler and CoRoT observa- 
tions from main-sequence (strong relation with T e $) to red-giant 
stars, there being a large dispersion in the temperature range 
4500±250K. Our results at lower temperature (r eff < 4200 K) 
again show that there is a strong dependence of the damping on 
r e ff. Damping measurements for stars in this range using for ex- 
ample the ongoing Kepler mission will demonstrate whether our 
description is valid for this range and put tighter constraints on 
the relation for red giants. 
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